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Noise and Interference 
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(Manuscript received September 13, 1968) 

A single, constant amplitude, in-band, additive interference is included 
in the analysis of detecting phase sh ift keyed signals in gaussian noise. For 
coherent detection we give a method applicable to any M-phase system, and 
evaluate the symbol error probability for M = 2, 3, and 4- For differential 
detection we treat the important cases M = 2, .'+, 8, and 16, offering com- 
prehensive numerical results for each. 

The analysis in each case is based on a single sinusoid with random 
phase adding to the noisy phase shift keyed signal. The results are then 
interpreted to include an angle modulation impressed on the continuous 
wave interferer. The receiver consists of an ideal phase discriminator with 
a perfect slicer. The channel is also assumed ideal in that intersymbol 
interference is not considered. 

I. INTRODUCTION 

Phase shift keying (psk) is becoming more popular as a modulation 
scheme for transmitting digital information. Lately much analysis has 
been done for both coherent and differentially coherent detection. 
Unfortunately the analyses done to date have generally considered 
only two signal degradations: channel anomalies (such as distortion, 
gain and delay variations, and so on) and thermally generated noise 
modeled by a gaussian random process. This article considers the 
effects of a spurious signal, or interference, falling in the band of the 
desired signal, as well as gaussian noise. It is understood that both 
the noise and the interference additively corrupt the desired signal; 
these are the only perturbing factors. 

For coherent detection, the phase probability density function for 
the received composite of signal, noise, and interference is found. 
From this, the theoretical error probability may be evaluated for 
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any iW-phase system. We give comprehensive numerical results for 
the important cases, M = 2, 3, 4. 

For differential detection we present analysis and results for M — 2, 
4, 8, and 16. In the binary case a simple closed form solution was 
found which yields both a good approximation and exact bounds to 
the actual error probability. The solutions to the multilevel M > 2 
differential detection problem, which are exact, required machine 
computation of a double integral; complete numerical results are given. 

Finally, we draw general comparisons between coherent and dif- 
ferential detection error performance as affected by interference. 

II. SIGNALS, NOISE, AND INTERFERENCE 

A phase shift keyed signal has the form (ignoring any amplitude 
function) 

s(t) = cos [2rf.t + 0.(0] (!) 

where we choose to normalize the peak signal amplitude to unity. 
The digital modulation is carried in the angle of s by <f> 8 {t), which 
assumes discrete values from a set of M equally spaced points in 
[0, 2tt] at the sample times T seconds apart. Thus the Nth. message 
or baud is modulated by 

*OT)-|f, * = 0,1,2, ••• ,M-1 (2) 

where each of the M values of k is equally probable. 

For a coherent receiver an M-ary symbol is transmitted in one 
baud by the value of k. For a differential detection receiver the in- 
formation is transmitted by the changes in k (or carrier phase) 
between adjacent bauds. 

The noise is presumed to originate thermally and is therefore 
modeled in the usual fashion by a stationary zero mean gaussian 
random process with uniform spectral density. At the output of a 
symmetrical bandpass filter the noise voltage may be written as 1 

n(t) = u(t) cos (27T/.0 - v(t) sin (2ir/,0 (3) 

where u and v are low-pass, stationary, independent, zero mean gaus- 
sian random processes with ensemble averages 

(n% v m (u\ v = (v 2 )„ = a 2 , (4) 

equal to the noise power. 
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In the differential detection analysis we make the further restric- 
tion that the noise process autocorrelation vanish at the baud inter- 
val, thus 

R nn (T) = 0. (5) 

This assures the four gaussian random variables 

u(t ), u(t + T), v(t ), v(t + T) 

to be uncorrelated and hence independent. 2 

Interference shall consist of a constant amplitude, possibly angle 
modulated, sinusoid which lies within the bandwidth of the detector. 
It is assumed to originate independently of the signal, and so it is 
natural that its phase relationship to the signal is random with all 
angles equally probable. Therefore let 

i(t) = b cos \2TTJit + <f>i(t) + T] (6) 

which has a peak value of b, and is angle modulated by <£ f . The arbi- 
trary phase angle r, independent of <j>i, is a random variable whose 
probability density function is (2-tt)- 1 when reduced modulo 2ir. 

For coherent detection, where the interference is observed only 
once per symbol, the random phase variable T vitiates the modulation 
<f>i because the sum (<f> t + V) is distributed exactly as if it were uni- 
form. This is discussed in Section III. 

III. COHERENT DETECTION 

An ideal phase discriminator is assumed which compares the re- 
ceived wave (composed of signal, noise, and interference) with the 
unmodulated signal carrier (the reference) and produces instantly 
the signed phase difference between the two inputs. 

The detector examines the discriminator output and announces 
an estimate of the transmitted symbol. The detector operates with 
no timing error and with zero width decision thresholds. Using maxi- 
mum likelihood detection based on equal a priori symbol probabilities, 
the thresholds are at ir/M, {Ztt)/M, . . . , [(2M - l)ir]/M. In a 
phasor diagram these thresholds correspond to (2tt)/M angular sec- 
tions centered about the M signal positions. 

The approach used to find Pe, the probability of a symbol error, 
is to find the probability density function of the phase of the re- 
ceived composite (s + n + i), and then integrate the density over 
the error regions. 
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We first notice that the phase angle of the interference relative 
to the signal is, from equations (1) and (6) , 

#(0 = 2tt(/,- - f,)t + *&) - <t>.(t) + r (7) 

where T is independent of the other terms on the right-hand side. 
Since r is uniformly distributed modulo 2rr, the relative interference 
phase process *(£) is also uniformly distributed modulo 2v. This is 
a general result for the modulo addition of several variables, one 
of which is uniform. 3 

Figure 1 is a phasor diagram of the receiver input components, 
signal, interference, and noise, at a sample time t . The phase refer- 
ence is 2irf s t + 0a so that the signal lies along the reference (vertical) 
axis. The orthogonal noise phasors are assumed to be at angles and 
it/2, relative to the signal. We seek the probability density function 
of the resultant angle A, and begin by considering the two dimen- 
sional joint probability density function of the cartesian coordinates 
of the resultant phasor. Conditioned on *, it is clearly jointly gaus- 
sian with means 

(x) av = b sin *, (y)„ = 1 + b cos $ (8) 

so that 

}xy(x, V I <t>) 

= 2^7 exp {" h [( * " b sin 0)2 + (y - 1 ~ bcos *)']}■ w 

Eliminating the <I> dependency gives 



exp {- ^ [.t 2 + (y - l) 2 + 6 2 ]} 



/*>■<*' tf- (2XCT) 



■f W exp fX [x 2 + (y- l) 2 ]* cos (0 + ,)} d<j> 

where rj = tan -1 [ (y — l)/x] is not a function of <f>. 
This integrates directly to 

hr{x > y) = 2b exp {~ h [x * + {y ~ 1)2 + 62] } 



(10) 



i \-2ix* + {y-i) 2 ]*} (ii) 
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-INTERFERENCE 




Fig. 1 — Phasor diagram of the signal, noise, and interference components at 
a sample time / . 

where 7 is the zero order modified Bessel function of the first kind. 
We now convert equation (18) to polar coordinates through the 
usual transformation 

x = r sin a and y = r cos a (12) 

which has the Jacobian r. Then the polar coordinate two-dimensional 
density is integrated over all radius values to yield the desired proba- 
bility density function of the angle. 

1a(oi) = 7^2 j exp < — ^5 [r + b 2 + 1 - 2r cos a]> 

■h\ \ (r 2 + 1 - 2r cos a) h \r dr. (13) 

The above integration has been done numerically to generate exact 
f.i(ct) curves for several values of 

-20 Iog,ol2 , '--'(T] = carrier to noise ratio in dB (CNR), 

—20 log 10 £> = carrier to interference ratio in dB (CIR). 

It is clear from equation (13) that f A (a) has at least the symmetries 
of cos a. 

Figure 2 offers typical families of f A probability density function 
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Fig. 2 — Probability density function of the phase of s + n + i for various 
CIR values, (a) CNR = 10 dB, (b) CNR = 15 dB, (c) CNR = 20 dB, (d) 
CNR = 25 dB. 
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curves for representative CNR values. As the interference amplitude 
increases it is seen to control or to affect the shape of the curve to 
a greater degree. At high interference levels a saddle-like shape ap- 
pears with peaks at roughly tan _1 b, as one would expect in the 
absence of noise. 

Since equation (13) is the probability density function of the 
angle of the complete receiver input, and since the probability of a 
symbol error is the probability that A lies outside the region \—ir/M, 
if/M] at t , we have 



Pe = [ '" /,,(«) da + f U(a) da 

J-T Jr/M 

which by symmetry is 



Pe = 2 



j A (a) da. 



(14) 



(15) 



Again, integral (15) was done numerically (a simple summation 
of the /,i data) for practical combinations of CNR, CIR, and for 
M = 2, 3, and 4. The results appear in Figs. 3-5.* 

The above method, which employs two (rather simple) machine 
integrations, yields exact results but at the expense of not generating 
useful expressions for Pe. Therefore we now indicate one approach 
which yields Pe for M = 2, and bounds Pe for M > 2, as a convergent 
series. We begin by considering binary reception. 

Referring back to Fig. 1, an error is made if | « | > 90° or, equiv- 
alently, if the resultant resides in the lower half plane y < 0. Then 
for fixed <J>, 



Pe | * = (• 



•' — z 



exp 



r-j (y - 1 - b cos 4>Y 

la 



dy 



= h erfc 



1 + b cos 

(2)V 



Averaging over the uniformly weighted * gives 



IT 



/ 1 + 6 cos \ . 
^ erfc, ~(2}V ,d0 " 



(16) 



(17) 



This integral, which is virtually the cumulative distribution func- 
tion of a sine wave of amplitude b plus gaussian noise of variance a 



* The abscissa values are true carrier-to-noise power ratios, and are not ad- 
justed to reconcile bandwidth to bit-rate differences. One may do this by sub- 
tracting 2 dB (3 dB) from the abscissa values for M = 3(4). 
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(evaluated at —1), has been examined by Rice 4 and others. It can be 
evaluated by expanding the integrand in a Taylor series about (2<r 2 ) -i 
and then integrating term by term. If the interference is small, b <C 1, 
only the first several terms need be retained for reasonable accuracy. 
The Pe values obtained for Binary may be used to bound the symbol 
Pe for M > 2. The decision thresholds are at ±ir/M for the ilf-ary re- 
ceiver. The error region consists of the union of two half planes formed 
by the extended detector thresholds. The probability that the resultant 
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Fig. 3 — Binary (M = 2) Pe versus CNR. Coherent detection. 
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Fig. 4 — Ternary (M = 3) Pe versus CNR. Coherent detection. 

phasor terminates in the error region, Pe, is thus bounded by the proba- 
bility of terminating in either half plane. As M increases, the bound 
becomes a good approximation* because the size (hence the relative 
probability) of the doubly counted intersection decreases rapidly with 
M. 6 By symmetry, the probability of terminating in either half plane is 
twice that of one half plane. 

* The approximation improves with CNR also. 
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Fig. 5 — Quaternary (M = 4) Pe versus CNR. Coherent detection. 

The probability of one half plane is related to the binary Pe very 
simply. The distance to the boundary from (0, 1) is sin (tt/M). If 
now the interference phasor and noise phasor were scaled by the 
same factor, we see that the probability of the half plane is just 
the binary Pe with interference amplitude b sin (w/M) and noise 
variance [<r sin (tt/M)] 2 * 



* The author is grateful to S. O. Rice for suggesting this notion. 
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Therefore we have the interesting relationship 

Pe (m > 2, a sin ^ , b sin ~\ £ 2 Pe (M = 2, a,-b) (18) 

which means that twice a given binary Pe value is an upper bound 
(or approximation) to an i\/-ary Pe where the binary CNR and CIR 
values are each increased by —20 log sin (n/M) dB. 

IV. DIFFERENTIAL DETECTION — BINARY 

This type of detection has become widely considered lately because 
it eliminates the requirement of phase synchronism between the 
transmitter and receiver. The price one pays for nonsynchronous 
detection, however, is poorer performance. 

The analysis for differential detection is complicated by the fact that 
the phase reference, being the previous signal, is subject to the same 
corruptions as the present signal being phase detected. (For the noise- 
only case, an exact solution in closed form is available for binary, and 
good approximations exist for M ^ 4.) We begin with an analysis for 
binary differential psk (d-psk), considering again a single CW inter- 
ference in addition to noise, which yields in closed form both bounds 
on Pe and a good approximation. 

We are concerned with the reception of two successive bauds, where 
the data is encoded in the phase change of the signal. Arbitrarily, let 
no phase change represent a "0" and a ■* phase change represent a 
"1." For convenience, we will refer to the signal during a baud 
interval as a "pulse" of carrier at a certain phase angle, although 
in a pure phase modulated (PM) system the signal would not con- 
sist of carrier pulses. 

From previous assumptions, the noise corrupting each signal pulse 
acts independently; the interference at two adjacent sample times 
U and U + T does not. This dependency may be summarized by an 
angle 

9= 2,r(/< -j.)T 

which is the relative phase slip of the interference from one sample 
instant to the next. Assume, for the present, that the interference 
modulation is absent. 

We will use the same equally spaced detection thresholds as in 
the coherent psk analysis. Ignoring interference, the probability of 
error in d-psk is not data dependent because of obvious symmetries. 
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However, the addition of interference destroys that symmetry, and 
causes Pe to be strongly data dependent for a given 6. Fortunately, 
however, the error probability for only one symbol (that is, "0" or 
"1") needs to be found because the probability of error for the other 
symbol (s) is derived directly from it. An over-all probability of 
error is then found by averaging the individual symbol error proba- 
bilities with equal weighting. 

Consider the transmission of a "0" whereby two carrier pulses of 
the same phase are sent. A "double exposure" phasor diagram, Fig. 
6, pictures the signal, noise, and interference components at the two 
successive sample instants t and U> + T. The interference at tn 
assumes an angle 4> relative to the signal, where £ is random and 
uniformly distributed in [0, 2tt] . At time U + T the interference has 
progressed to an angle <f> + 9. 

The noise phasor amplitudes are the random variables 

U d m u(t ), V d = v(t ), U = u(t + T), and V m v(t + T) (19) 

which we recall are independent, equal variance, zero mean gaussians. 
The two resultant phasors, Z and Z d , are the actual phase dis- 
criminator inputs; the output being their phase difference, 8. Since a 



INTERFERENCE 



INTERFERENCE 




Fig. 6 — Signal, interference, and noise phasors at successive sample times 
to and U + T for a transmitted "0." 
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"0" is transmitted, the probability of error is 



Pe I "0" = Pi 



>I 



— 7T < 5 ^ 7T 



which is equivalent to Pr{cos 8 < 0}. 
Also, because 



cos 5 = 



Re [ZZ%] 



(20) 



(21) 



,z\\z d \ 

and the denominator is nonnegative, we have alternatively 

Pe | "0" = Pr{Re [ZZ* d ] < 0}. (22) 

We now employ a technique used by Stein which leads to a closed 
form solution to equation (22) . 6 Based on the simple identity 



(23) 



z + z d 


2 


z - z d 


2 




2 



Re [ZZ* d ] = 
we see that 

Pe| "0" = Pr{|Z + Z d \ < \Z - Z d \) 
For economy of notation we let 

2 = Z + Z d = U -f U d + i(V + V d ) + C(S) 

A ■ Z - Z d = U - U d + i(V - V d ) + C(A) 
where C( ) denotes the nonrandom components,* that is, 

C(S) = 6[sin + sin (<f> + 6)] 



(24) 



(25) 



(26) 



+ i{2 + 6[cos tf> 4- cos (<f> + &)]} 
C(A) = 6[sin (0 + &) - sin <£] + ib[cos (<f> + 6) - cos 0] . 

Thus 2 and A are complex (two dimensional) random variables 
whose jointly gaussian orthogonal components have equal variance 
2<t 2 . We are concerned with the magnitudes, | 2 | and | A |, often 
referred to as a Ricean random variable. 

The probability density function (of | 2 |, for example) is the well 
known result for the envelope of sine wave plus noise 4 



4 m r r ^4- ices) i 2 "ur nc(s) n 



(27) 



* is constant and we have conditioned the solution on tf>. 
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Furthermore, 1 5 1 and | A | are independent by the following 
argument. U + U d is independent of V — V d . Consider U + U d and 
U — U d . They are uncorrected, hence independent, by virtue of 
their equal variance. 

<[(£/ + U d )(U - U d )])„ = [U 2 - U 2 d ] = (28) 

Therefore all four random components of 2 and A are independent, 
from which the independence of | 2 | and | A | follows. 

The advantage of this approach is that equation (24) , the proba- 
bility that the amplitude of one complex gaussian (that is, a Ricean) 
exceeds another, is expressible in terms of the function 7 * 

Q(A, B) = f r exp (- — + T ' )/ (Ar) dr. (29) 



One formulation is 8 



means from equation (2( 
C(2) | = 2J~1 + 26 cos | cos (* + f ) + & cos 2 |J 



Evaluating the magnitudes of the means from equation (26), 

11 



(30) 



(31) 



C(A) | = 2&sin^- 



We have thus far found Pe|"0" exactly, conditioned on the initial 
interference angle <f>. The desired result is obtained by averaging 
equation (30) over <£. However, inserting equation (31) into equation 
(30) leaves an expression which offers little promise for analytically 
performing the integration. As an alternative, the integral is both 
approximated and bounded in what follows, thereby avoiding a ma- 
chine integration. 

The integration parameter <£ appears in |C(2)| but is not a factor 
of I C(A) I. Then with the identification 



B - | C(2) \/2a 
A = | C(A) |/2<r 



(32) 



* Q functions are tabulated, but not adequately in the argument ranges needed 
for these problems. Their usefulness lies in having good approximations which 
lead to easy machine calculations. See Ref. 8. 
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we are led to consider the behavior of the integrand 

F(A, B) = Ml - Q(B, A) + Q(A, B)} (33) 

as a function of B only (fix A). Familiarity with F(A, B) suggests 
that it is roughly exponential in —B-/2. 
To affirm this supposition we first show that 

dF( ^ B) = ~BK(A , B)F(A , B) (34) 

where K(A, B) is bounded near 1: 

1 ^ K(A, B) ^ 1 - e. (35) 

This approach is motivated by the recognition that if F(A, B) 
could be approximated by an exponential in B 2 , the integral of 
F(A, B) over <f> would be simple. 

The derivation of equations (34) and (35) consists of elementary 
manipulations of the series and integral representations for Q(A, B) 
which are given in Ref. 8. Henceforth, a prime designates partial 
differentiation with respect to B. For brevity, let 

£^exp[- A '+*' ]- (36) 

Then using the series representations 

Q{A,B) =EJ:(^j l I m {AB) (37) 

1 - Q(B, A) =EZ (£fl m (AB) (38) 

we see that 

2F(A, B) = 2Q(A, B) - EI (A, B). (39) 

Then 

2F'(A,B) = 2Q'(A, B) - EI' (AB) - E'I (AB) (40) 

= 2Q'(A, B) - EAI^AB) + BEI (AB). (41) 

Referring to the integral representation for the Q function we 
inspect 

Q'(A, B) = -BEh(AB). (42) 
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Therefore, from equations (42) and (41) we have 

2F'{A,B) = -BE^UAB) +|/ 1 (A5)]- 



(43) 



But from equations (37) and (38) , 



1 -Q(B,A)+Q(A,B) =E 

So that equation (34) is true for 



I Q {AB) + 2 £, (|) m / m (^B)] • (44) 



I (ABj + (4Vi(^) 

K(A, B) = — XW - (45) 

I Q {AB) + 2[^)I 1 (AB) + 2 £ (^J Z m (AB) 

Since .4, £, and the I m (x) are nonnegative the unity upper bound 
in equation (35) is obvious. For the lower bound, invert equation 
(45) and notice that 

S 1 + 2 t (I)" (47) 

since the 7 m (re) decrease with the order m. 

Referring to the definitions of equations (32) and (31) it can be 
shown that A < B whenever b < l/(2) % , so that the summation con- 
verges for reasonable interference levels (CIR > 3 dB). Then 

K{A,BT l £1+^2 (48) 



so that 

- 2A 



B + A 



(49) 



In most cases B ~5> A so that K(A, B) is near 1. 

We now solve the linear, first order differential equation (34) to 
obtain 



F(A,B) = Pe exp 



[ - tK(A, t) dr 



(50) 



with initial data F{A, B ) = Pe . The above expression for F(A, B) 
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is exact. We now approximate the exp argument by approximating 
K(A, B) with a constant K . Then we are able to carry out the 
integration in equation (50) to obtain 



F(A,B) ^Pe exp 



^(B 2 - />' 



»]• 



(51) 



We know the function at B , and we extrapolate to the function 
using an approximation of its derivative. Now the <f> integration will 
range over some section of the exponentially varying F{A, B) in 
equation (51) above. Because it is exponential, the significant con- 
tribution to the integral occurs over a relatively small range of B 
where F(A, B) is near its maximum. This suggests that the initial 
data be specified at a maximum so that the approximation function 
is best in the important range of integration. Therefore we let 



B = J5 miI 



so that 



Pe = F(A, B ) ^ F(A, B), 
We notice from equation (31) that 



< <b < 2tt. 



SmlB ™ - ( 1 — 6 COS ^ 



and so 



n2 „ 2 26 6 
B 2 - Bl = -a cos ■= 



1 + cos ( 4> -f- 2 



(52) 



(53) 



(54) 



(55) 



In addition, we slope-match at B , approximating K(A, B) by K = 
K{A, Bq). From equation (45) 



1 + 



A\ UABo) 



Bo) h(ABo) 



. A^KABo) 



K(A,B Q ) 

1 + 2 \BJ I (AB ) 
Inserting equation (55) into equation (51) and integrating, 

£ f T F(A, B) d<f> - F(A, Bo) exp [- ^ cos |] 

i 



(56) 



2^ eXp 



f^cos^ + l 



(57) 
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gives the desired result 
Pe | "0" = i[l - Q(B , A) + Q(A, B )] 



•exp 



f bK e~\ T [bK n e~\ . rQ , 

L"^" cos 2j /o L^ cos 2_r (58) 



Exact bounds are now easily obtained by bounding K(A, B). With 
B still chosen to be B min , clearly 

[ - tK(A, t) dr ^ f - rK min dr = -\{B 2 - Bl)K mia (59) 

where /v mln is the minimum value of K(A, B) in [B , B]. Similarly, 



Therefore 



so that 



f - tK(A, t) dr ^ -\{B 2 - Bl)K max . 



F(A, B) -^ Pe exp [-h(B 2 - Bl)K max ] 



(GO) 



(61) 



[-K^Aeosl] 



IN | ()", 0-^-Pe o exp | -i^ min -j cos ^ |/ \K min \ cos | ■ 



a 2 cos - | • (62) 



The ratio of the bounds in equation (62) may be bounded in order 
to ascertain their closeness. We omit the cumbersome derivation, but 
state the result below. 



Upper Bound 
Lower Bound 



< 



1 0,0 

- - cos - + sin - 



7 — cos - — sm - 
oil 



(63) 



For CIR ^ 10 (15) dB, the ratio is less than 1.6 (1.25). 

The symbol "1" is transmitted as two pulses of carrier 180° apart. 
Then correct reception results if 90° < | 8 \ < 180°. If an analysis 
quite similar to the preceding were carried out for this case, the 
resulting expressions would be identical with those above except that 
is replaced by it — 0. It follows that 

Pe | "1", = Pe | "0", it - 0. (64) 

We will elaborate on the relationship of Pe for the different data 
symbols in Section V. 
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Notice that Pe|"0" is symmetric in 6 about and v. The sym- 
metry follows from the averaging of <f>, and the insignificance of 
which of the two pulses occurs first as far as the detector is concerned. 

Since the overall error probability 

Pe (6) = i[Pc | "1", 6 + Pe | "0", 6] (65) 

we use equation (64) to write 

Pe (6) = |[Pe | "0", + Pe | "0", r - 0] (66) 

which is easily shown to be evenly symmetric about 0, x and tt/2. We 
therefore need examine Pe only in the range 0^0^ 90°. 

Pe(0) does vary considerably as seen in Fig. 7. Here the maximum 
and minimum values of Pe(0), which happen to occur at -tt/2 and 
respectively, are given for interesting combinations of interference 
and noise levels. To further illustrate the effects of interference, we 
present curves of decibel degradation versus in Fig. 8. Degradation 
is defined as the dB reduction in carrier-to-noise ratio which is 
allowed to maintain the same Pe after removing the interference. 

Finally, consider an angle modulation impressed on the interfer- 
ence. This situation may be viewed simply as a time varying 0. Then 
one may average the Pe(0) results given here over the variations of 0. 
If this is undesirable, the curves of Fig. 7 are certainly bounds on 
Pe averaged over the variation. 

V. DIFFERENTIAL DETECTION QUATERNARY 

We now examine the effect of a single interference on a differentially 
detected quaternary (4-phase) signal. We will refer to the four sym- 
bols as "0", "1", "2", and "3", where the associated baud to baud 
phase shifts are 0, tt/2, it, and —tt/2 respectively. As before, the 
phase discriminator examines the two composite phasors, Z and Z d , 
and reports their angle difference 8. The ordering of the bauds is 
important, since we must distinguish between 8 = tt/2 and 8 = — n/2. 
The "0" symbol possesses the same symmetry as in the binary case; 
we therefore base the analysis on Pe|"0". Then in an analogous fashion 
we relate Pe|"0" to the probability of error for the other symbols. 

In the binary case the receiver tested for the sign of Re [ZZ*,]. This 
test was transformable to a test between the amplitudes of two Ricean 
random variables, one which enjoys a closed form solution. Unfortu- 
nately, in the present case the test, which is for a "0" 
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Pe | "0' 



. Re [ZZ*] 
Pr<cos5 = 



= 2 



is nob known to be transformable to one which has a closed solution. 
On the other hand, we offer a very straightforward analysis which 
is exact and amenable to machine computation. 

Figure 9 is a "double exposure" of the signal, noise, and interfer- 
ence components for a "0". The two carrier pulses are, of course, coin- 
cident and lie along the reference axis. We recall that the angle of 




10 12 14 16 

CARRIER -TO- NOISE RATIO IN DECIBELS 



Fig. 7 — Binary Pe versus CNR. Differential detection. 
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the interference at t is <£, and at t + T is <f> + 6. We will again 
average over </>, leaving 6 as a parameter. 

At time t the resultant of signal and interference, shown by a 
dashed line, has length W and angle £ which by inspection are 



W = (1 + b 2 + 26 cos 0)* 

-i T b sin ^ 1 

£ = tan iir. * 
LI + b cos <£j 



(G7) 
(68) 




8dB 

|_ CARRIER-TO- 
NTERFERENCE 
RATIO 



10 20 30 40 50 60 70 80 90 
IN DEGREES 



Fig. 8 — Degradation in CNR caused by interference versus for binary dif- 
ferential detection. 
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NOISE 



-NOISE + INTERFERENCE 
CORRUPTED FIRST PULSE 




Fig. 9 — Phasor diagram for quaternary differential detection. 

To the resultant of signal and interference adds a Rayleigh ampli- 
tude, uniform angle noise phasor. The resultant W is perturbed both 
in amplitude and angle by the noise. The resulting amplitude is 
unimportant, but the angle, | + a, establishes the reference for de- 
tecting the second pulse. The probability density function of <x is well 
known to be 



/«(«) = 7T ex P ( — *) I * + *** cos a 

-IT 



•exp (^ cos 2 a)[l + erf (** cos a)] 



where 



2 r 1 

erf (x) s -r / exp (-w 2 ) 

7T Jo 

is the usual error function integral and 

W 2 



da 



y = 



2<x 2 



(69) 



(70) 



(71) 



is the power ratio (signal plus interference) /noise. 
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At t + T the second pulse is examined. It is disturbed by inter- 
ference to the point (.r , ijo) , 

.!•„ = b sin (0 + 6) , 7 2) 

y = 1 + b cosfa + 8). 

The decoding region is a quarter plane bisected by the direction 
of the first pulse as shown by the orthogonal decision boundaries. 
The probability of correct reception is simply the probability that 
a random noise phasor originating at {xo, y ) will terminate inside 
this quadrant. Using the independent orthogonal phasor representa- 
tion for noise, and choosing the components to lie alongside the per- 
pendicular distances di and d 2 from (x , y ) to the boundaries, we 
write directly 

1 - Pe | "0", *, « - | [l + erf (gfc)] \ [l + erf gfc)]. (73) 
Now using erf + erfc = 1 we have 

Pe | "0",4>,a = |erfc(|y 

+ , P1 ,,(^)_,erfc(^)erfc(^). (74) 
The distances <U and </•_. may be verified to be 

d, = -y sin I a + £ - | J + a-'o cos (a + $ - |J 

d 2 = ?/„ sin ^* + £ + |j - x a cos ^a + 5 + |J 

such that they take the positive sign if (.r , y ) lies on the correct 
reception side of the respective boundary. 

Eliminating the <f> and a dependency results in a finite limits double 
integral 

Pe I ••0" = f ^ f 1M Pe | "0", 0, a da (76) 

' Jo 2* J-T 

which was machine evaluated. 

The relationship between Pe | "0" and the other symbols is easily 
demonstrated graphically. Figure 10 is a phasor diagram illustrating 
a typical noise and interference corrupted first pulse, and the four 



(75) 
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ARBITRARY DIRECTION 

OF CORRUPTED 

FIRST PULSE 




Fig. 10 — The relationship of Pe\0 for the individual symbols. 

possible positions of the pre-interference second pulse. Decoding 
quadrants determined by the angle of the first pulse are shown by 
dashed lines. Adding to each of the possible second pulse carrier 
phasors, which correspond to the four data symbols, is an interfer- 
ence phasor having some angle <f> + 8. 

If a "0" had been sent, the probability of error conditioned on the 
present geometry (that is, the values of <£, a, and 6) is the probability 
that a noise phasor originating at the tip of the solid interference 
phasor (labeled "0") terminates outside the "0" quadrant. This proba- 
bility is a function of a and the distances from the tip of the inter- 
ference phasor to the boundaries. 

Now assume that a "3" had been sent instead. The associated 
second carrier pulse is shifted clockwise by it/2, but the interference 
phasor at ta + T is still at <f> + 9 relative to the diagram reference. 
Returning to the "0" quadrant, consider an interference phasor hav- 
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ing the angle <f> + 6 — -rr/2, shown by the dashed line labelled "3". 
It is apparent (and trivial to show) that the distances of the dashed 
"3" phasor to the "0" quadrant boundaries are identical with the 
distances of the actual "3" interference phasor to the "3" quadrant 
boundaries. We conclude, therefore, that 

Pe| "3",0,a, = Pe| "()»,*,«, 0-|- (77) 

Now integrating both sides of the equality over all <f>, a yields the 
desired relationship 



Pe | "3", = Pe | "0", - f • 


(78) 


Similarly we have 




Pe | "2", = Pe | "0", 6 - t 


(79) 


Pe | "1", = Pe | "0", + | 


(80) 


so that the average symbol error probability becomes 




Pe (0) = ! Pe | "0", + Pe | "0", - | 





IV "(>". o + I + p e | "0", - 7T (81) 



2 



which is solely in terms of Pe | "0". 

Notice that equation (79) is exactly the result obtained for a "1" 
in binary. This is not surprising, since a "2" constitutes a 180° phase 
shift of the second pulse. In fact, the arguments relating to Fig. 10 
may be generalized for an il/-phase d-psk signal with a "J" symbol 
phase shift of (2irJ)/M, to be 

Pe | "J", = Pe | "0", d-\- l ~- (82) 

Averaging over the symbols in equation (81) produces a Pe(0) 
which is evenly symmetric about 0, tt/4, tt/2, . . . ; points half as far 
apart as in binary. Also, while Pe|"0" still varies over a considerable 
range, when four symbols are averaged instead of two the range of 
Pe(0) is significantly decreased. This is evidenced by the numerical 
results plotted in Fig. 11. Again, the solid and dashed lines represent 
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Fig. 11 — Quaternary symbol Pe versus CNR. Differential detection. 



maximum and minimum values of the symbol Pe. The values of 6 
which correspond to the maximum and minimum are not the same for 
all noise and interference level combinations. However, Pe(0) is gen- 
erally lowest near 0, r/2, it . . . and highest near ir/ 4 > 37r / 4 > • • • » 
and so on. 

Since Pe (8) fluctuates less severely for quaternary, it is more mean- 
ingful to average over 6. This was done, and the average used as a 
base for computing the degradation curves of Fig. 12. 
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VI. DIFFERENTIAL DETECTION — M > 4 

For large M systems, for example M ^ 8, we compute the average 
error probability 



Pe = — / Pe (0) dd 

7T Jo 



(83) 



assuming now that is a uniformly distributed random variable. Be- 
cause of the rapidly diminishing 8 dependency noted in Section V, 
the average over 6 is a useful measure of error performance. The 
choice of a uniform distribution for 6 allows an approach which relies 
on the previously obtained /a(«) data, rather than finding Pe(0) and 
then integrating. 

Let the interference phase angles at t and / + T be random variables 
$ and<3> -f- 0, respectively. We first note that the sum <S> -f- modulo 2t 
is uniformly distributed since $ (or 0) is uniform. Hence both inter- 
ference angles, $ and $ + 0, are uniformly distributed. Furthermore, 
since for any $ = <£ the sum <f> + (modulo 2if) is uniform, we conclude 
that «£ and <f> + are independent. The adjacent interference angles 
then are independent; and consequently the phase angle, A, of s + n + i 
is independent from sample to sample. 

We therefore obtain the probability density function of the differ- 
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Fig. 12 — Degradation in CNR caused by interference for quaternary differen- 
tial detection. 
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ence angle 8 = A{t + T) - A(t ) as the modulo 2tt convolution of 
j A {a) with itself. Then the integral of the probability density func- 
tion of 8 over \S\ > (v/M) yields Pe. Numerical results were obtained 
in this fashion for M = 4, 8, and 16. The M = 4 data was in excellent 
agreement with Fig. 11. Pe for M = 8 and M = 16 is displayed in Fig. 
13. When b > sin(ir/2Af) the interference alone can exceed the 
thresholds and cause errors. This is seen as a Pe floor for low CIR 
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Fig. 13 — Symbol Pe versus CNR for 8 and 16 phase differential detection. 



PSK ERROR PERFORMANCE 441 

values; increasing CNR does not cause Pe to tend toward zero. 

Finally, we remark that the averaging of Pe over the M symbols 
was implicitly done in the 6 averaging. We see from equation (82) 
that all symbols have equal error probability when 6 is uniformly 
distributed. 

VII. SUMMARY AND CONCLUSIONS 

We have evaluated the symbol error probabilities for both coherent 
and differential detection of low M psk signals in the presence of in- 
terference. The effect of the interference is readily observed in the 
curves of Figs. 3-5, 7, 11, and 13. Although the increase in Pe result- 
ing from interference is large, it is considerably less than if the inter- 
ference were replaced by gaussian noise of the same power, especially 
at low Pe levels. 

Comparing Figs. 3, 4, and 5, we see that vulnerability increases 
with M. That is, at a given CNR, the Pe is raised least for M = 2 
and greatest for M = 4 by the addition of interference. For example, 
without interference the error performance of M = 2 versus M = 4 
differs by 3 dB. When a —10 dB interferer is added, it differs by ap- 
proximately 5 dB, indicating a 2 dB CNR penalty for equal Pe values. 

Drawing comparisons between coherent and differential detection 
reveals that differential detection clearly suffers more degradation. 
Binary differential, however, performs about as well with interference 
at optimum 8 values as does binary coherent. This is in contrast with 
the performance disparity between the two M = 4 systems. With a 
— 10 dB interference, differential detection suffers a degradation rang- 
ing from 5% to 8 dB; coherent detection is degraded only 4 to 4% 
dB for the same CNR range. 

We use degradation rather than the raw Pe versus CNR curves to 
make the above comparisons because of the inherent difference in 
performance between differential and coherent psk for noise alone. 
That is, the degradation comparisons automatically reconcile any 
disparities in the noise-only performances of the various systems. 

REFERENCES 

1. Davenport, W. B. nnd Root, W. L., Random Signals and Noise, New York: 

McGraw-Hill, 1958, pp. 158-159. 

2. Papoulis, A., Probability , Random Variables and Stochastic Processes, New 

York: McGraw-Hill, 1965, pp. 182-183. 

3. Scire, F. J., "A Probability Density Function Theorem for the Modulo y 

Values of the Sum of Two Statistically Independent Processes," Proc. 
IEEE, 56, No. 2 (February 1968), pp. 204-205. 



442 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1969 

4. Rice, S. 0., "Statistical Properties of a Sine Wave Plus Random Noise," 

B.S.T.J., 27, No. 1 (January 1948), pp. 109-114. 

5. Arthurs, E. and Dym, H., "On the Optimum Detection of Digital Signals in 

the Presence of White Gaussian Noise — a Geometric Interpretation and 
a Study of Three Basic Data Transmission Systems," IRE Trans, on 
Communication Systems, CS-10, No. 4 (December 1962), pp. 336-372. 

6. Stein, S., "Unified Analysis of Certain Coherent and Noncoherent Binary 

Communications Systems," IEEE Trans, on Information Theory, 1T-10, 
No. 1 (January 1964), pp. 43-51. 

7. Marcum, J. I., "Tables of the Q Function," Rand Corporation Memorandum 

RM-399 (January 1950). 

8. Schwartz, M., Bennett, W. R., and Stein, S., Communication Systems and 

Techniques, New York : McGraw-Hill, 1966, Appendix A. 



